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Measurement of Z boson Production in Lead-Lead Collisions at 



The ATLAS experiment has observed 1995 Z boson candidates in 0.15 nb~' of integrated luminosity obtained in the 201 1 LHC 
Pb+Pb run at ^s NN =2.16 TeV. The Z bosons are reconstructed via di-electron and di-muon decay channels. The results from 
both channels are consistent with each other and are combined. The background is less than 3% of the selected sample. Within the 
statistical and systematic uncertainties, the per-event Z boson yield integrated over rapidity \y z \ < 2.5 is proportional to the number 
of binary collisions estimated by the Glauber model. The elliptic flow coefficient of the azimuthal distribution of the Z boson with 
respect to the event plane is consistent with zero. 
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1. Introduction 

Extensive heavy ion (HI) programs carried out by the experiments at the Relativistic Heavy Ion Collider (RHIC) 
at BNL, and the Large Hadron Collider (LHC) at CERN, have established that the hot and dense matter produced in 
high energy collisions of HI imposes a significant energy loss on the energetic color charge carriers penetrating such 
a medium E[2)- An understanding of this phenomenon requires measuring the unmodified production rates of the 
particles in HI collisions before they lose energy. The best candidates to perform such measurements are particles 
which do not interact via the strong force. 



This note presents a measurement of Z boson production in Pb+Pb collisions at y/Sfw = 2.76 TeV using the 
electron and muon decay channels. The results demonstrate that the Z boson production rate in HI collisions is 
proportional to the mean number of binary nucleon-nucleon collisions ((N co n)) calculated within the framework of 
the Glauber model [3|. The analysis of the elliptic flow coefficient V2 of the Z boson angular distribution reveals no 
correlation between the emission angle of the particle and the azimuthal direction of the event plane. 

The ATLAS detector [4] at the LHC covers nearly the entire solid angle around the collision point. It consists of 
an inner tracking detector surrounded by a thin superconducting solenoid, electromagnetic and hadronic calorimeters, 
and a muon spectrometer incorporating three superconducting toroid magnets. 

The inner-detector (ID) is immersed in a 2 T axial magnetic field and provides charged particle tracking in the 
pseudorapidity range I77I < 2.5. A high-granularity silicon pixel detector covers the vertex region and it is followed by 
a silicon microstrip tracker and a transition radiation tracker. 

The calorimeter system covers the range I77I < 4.9. Within the region I77I < 3.2, electromagnetic calorimetry is 
provided by barrel and end-cap high-granularity lead liquid-argon (LAr) calorimeters, with an additional thin LAr 
presampler covering \r/\ < 1.8. Forward calorimeters (FCal) are located in the range 3.1 < I77I < 4.9. 

The muon spectrometer (MS) comprises separate trigger and precision tracking chambers measuring the deflection 
of muons in the magnetic field generated by superconducting air-core toroids. The precision chamber system covers 
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the region \rj\ < 2.7 with three layers of monitored drift tubes (MDT), complemented by cathode strip chambers (CSC) 
in the innermost layer of the forward region, where the background is highest. The muon trigger system covers the 
range \rj\ < 2.4 with resistive plate chambers in the barrel, and thin gap chambers in the end-cap regions. 



2. Analysis 

This analysis uses 201 1 LHC lead-lead collision data collected with the ATLAS experiment and corresponding to 
an integrated luminosity of approximately 0.15nb~'. 

Electron candidates were identified at the first trigger level (LI) as a cluster formed with (Ac/> x A77) = 0.1 x 0.1 
trigger towers of the electromagnetic calorimeter, covering the pseudorapidity range I77I < 2.5, excluding the transition 
region between the barrel and end-cap calorimeters (1.37 < I77I < 1.52). Candidates were then identified using the 
standard ATLAS reconstruction algorithm 0, requiring the matching of a track to an electromagnetic energy cluster. 
In addition to the track matching requirement, cuts based on the balance between the track momentum and cluster 
energy (E/p) and shower shape variables, were also used for electron identification. Electron selection is limited to 
I77] < 2.5 and both electrons are required to have E T > 20GeV. 

Muon candidates were selected using all three levels of the trigger system. The LI muon trigger searches for 
patterns of hits consistent with muons of a certain pj in the trigger chambers within I77I < 2.4. In addition, a full scan 
muon reconstruction was performed by the high level trigger to identify muons with pj > lOGeV. 

Single muons were reconstructed with varying levels of quality |6|. High quality muons were reconstructed in 
both the MS and ID subsystems with consistent angular measurements, as well as with a good match to the event 
vertex. At least one muon in each pair, matched to the trigger, is required to be of such quality. If the second muon in 
the pair has hit patterns in the MS and ID subsystems satisfying criteria of high reconstruction quality, the minimum 
pj threshold is set to lOGeV on both muons. If the second muon fails this condition, both muons were required to 
satisfy p T > 20 GeV. 

Once identified, electron and muon candidates were then paired in ee or pp. opposite charge combinations. Pairs 
with an invariant mass between 66 and 102 GeV were accepted as Z boson candidates. The Z boson sample has a 
contamination from background processes of approximately 5% in the Z — > e + e~ channel and 1% in the Z — > p + p~~ 
channel, mostly from combinatorial sources. The invariant mass distributions of the selected events together with 
estimated combinatorial backgrounds are shown in Fig. [T] Backgrounds from electroweak processes and top pair 
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Figure 1. The invariant mass distributions of Z — •> e + e~ (left) and Z — » (right) in data (points) and simulation (histogram) integrated 
over momentum, rapidity, and centrality. The simulation is reweighted to match the centrality distribution in data and normalized in the region 
66 GeV < mu < 102 GeV (/ = e,/i) indicated by dashed lines. Numbers of counts given in the plot correspond to the same mass region. 



decays JT] are small compared to the combinatorial background, and their contribution is neglected. The Z recon- 
struction efficiency is based on study of Pythia (version 6.425) [8] Z — > Z + T events with 66 GeV < m z < 116GeV 
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and \y z \ < 2.5 embedded into Pb+Pb events generated by the Hijing (version 1.38b) [9] event generator. The response 
of the ATLAS detector to the generated particles was modeled using GEANT4 irTOlfTTI . 

The main sources of systematic uncertainties on the measured yields in both lepton channels are associated 
with the precision to which the corrections applied to the data can be calculated. The efficiency and resolution 
uncertainties depend on the Z boson kinematics and are on average 8% (5.5%) and 2.5% (2.5%), respectively for 
Z — > e + e~(Z — > yu + /i~). 

Another way to search for any effect on Z — > l*l~ production due to the medium is to study the elliptic flow, v^> lfT2l . 
of the Z boson. The V2 is defined as the amplitude of the second Fourier harmonic of the Z boson azimuthal emission 
angle with respect to the event plane. The event plane contains the momentum vectors of the colliding nuclei, and is 
measured on an event-by-event basis from the azimuthal distribution of energy deposition in the FCal detectors. 



3. Results 

For both Z — > e + e~ and Z — > fi + p^ analyses, correction factors to account for the efficiency and detector resolution 
within the selected acceptance based on the simulation are calculated differentially in event centrality, p z , and y 2 . In 
each decay channel, the correction factor is applied and the corrected background, estimated by the same-sign pairs, is 
subtracted. The two decay channel measurements are then combined with weights set by their respective uncertainties. 
The fully corrected p 2 , and y z distributions for the Z boson are shown in Fig. [2] The data distributions plotted in 
(0-80)% centrality agree in shape with the Pythia simulations of Z boson production in pp collisions. 
The Z boson yields per-event, divided by (iV co ii), are shown in Fig. [3] The figure demonstrates that the Z — > e + e~ and 
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Figure 2. CoiTected transverse momentum (left) and rapidity (right) distributions of the measured Z — > I + l (I = e,/i) data compared to Pythia 
predictions normalized by area. The lower panels display the ratio of data over Pythia. 



Z — > p + p~ results are consistent as a function of (N part ) within their uncertainties for all p z and centrality regions. A 
combined uncertainty, including common uncertainties on the determination of (N co n) is presented with the horizontal 
brackets. Within the statistical significance of the data sample, the Z boson per-event yield obeys binary scaling. The 
dashed lines represent constant fits to the combined data as a function of the mean number of nucleons participating 
in the event, (A^ part ). 

The v 2 of the Z boson in the (0-60)% centrality interval is measured to be -0.015 + 0.018(stat.)+0.014(sys.). This 
observation is an independent measurement consistent with Z — > l + V yields being unaffected by the medium in HI 
collisions. 
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Figure 3. Centrality dependence of Z boson yields divided by (N c0 \{), measured in |v z | < 2.5. Results for ee and fifi channels are shifted left and 
right respectively from their weighted average, which is plotted at the nominal (N^ mt ) value. The statistical (bars) and systematic (shaded bands) 
uncertainties are calculated using the appropriately weighted average of the two contributing sources. Brackets show the combined uncertainty 
including the uncertainty on (N co \\). The dashed lines are constant fits to the combined results. 



4. Conclusions 

Z boson production has been measured in Pb+Pb collisions at -^1^=2.76 TeV using 0.15 nb -1 of integrated 
luminosity collected in the 2011 LHC run. Within \y z \ < 2.5, and 66 GeV < m u < 102 GeV, a total of 772 and 
1223 Z boson candidates are reconstructed in the Z — > e + e~ and Z — > jif + /i~ channels respectively. The combinatorial 
background is at the level of 5% in the electron channel and 1% for the muon channel. Yields of the Z boson 
production integrated within |y z | < 2.5 are consistent in the two channels in all measured pi and centralities. The 
momentum and rapidity distributions of the Z bosons measured in (0-80)% centrality region are consistent in shape 
with the Pythia simulations of Z boson production in pp collisions. Within the uncertainties the Z boson yield is 
found to be proportional to (N co \i). The elliptic flow coefficient of the Z boson is consistent with zero within the 
uncertainties of the measurements. 
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